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Plasma Contactors for Use with Electodynamic Tethers for Power Generation 

D.E. Hastings 1 and N. A. Gatsonis 2 
Dept of Aeronautics and Astronautics,MIT, Cambridge, MA 02139 

Abstract 

Plasma contactors have been proposed as a means of making good electrical contact between 
biased surfaces such as found at the ends of an electrodynamic tether and the space environment. 
A plasma contactor is a plasma source which emits a plasma cloud which facilitates the electrical 
connnection. The physics of this plasma cloud is investigated for contactors used as electron 
collectors. The central question that is addressed is whether the electrons that are collected by 
a plasma contactor come from the far field or by ionization of local neutral gas. This question is 
important to answer because the system implications are quite different for the two mechanisms. It 
is shown that contactor clouds in space will consist of a spherical core possibly containing a shock 
wave. Outside of the core the cloud will expand anisotropically across the magnetic field leading 
to a turbulent cigar shape structure along the field. This outer region is itself divided into two 
regions by the ion response to the electric field. A two dimensional theory of for the outer regions 
of the cloud is developed. The current voltage characteristic of an Argon plasma contactor cloud 
is estimated for several ion currents in the range of 1-100 Amperes. It is suggested that the major 
source of collected electrons comes by ionization of neutral gas while collection of electrons from 
the far field is relatively small. 
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Magnetic field 
ion acoustic velocity 
electron diffusion coefficient 
ion diffusion coefficient 
Electric field measured in frame moving with source 

Motional electric field 
ionization energy 
initial ionization fraction 
total current= I e + I c 
inflowing electron current 
outgoing ion current 
density scale length 
•>: 4 T- electron mass 

ion mass 
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f ambient ion density 

electron density 
contactor ion density 
neutral density 
radial distance 
radial dimension of source 
radial dimension of outer boundary of transition region 

radial dimension of core 
electron temperature in energy units 
ion temperature in energy units 
neutral, temperature in energy units 
Velocity of plasma source 
Alfven velocity 
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V D 

Differential electron ion drift velocity 

V n 

neutral velocity 

Vfhe 

Electron thermal velocity 

Vthi 

Ion thermal velocity 

fi 

Plasma diamagnetic parameter 

K 

electron thermal conductivity 

Vt 

Total electron momentum exchange frequency 

Vei 

Electron ion momentum exchange frequency 

Vtn 

Electron neutral momentum exchange frequency 

^ acoustic 

Electron momentum exchange frequency due to ion acoustic turbulence 

^ Bunt man 

Electron momentum exchange frequency due to Buneman turbulence 

<t> 

Potential = -V<£) 

<f>0 

Potential at center of cloud(r = ro) 

Up* 

electron plasma frequency 

n e 

Electron gyrofrequency 

n < 

Ion gyrofrequency 

Pe 

Thermal electron gyroradius(= t >the/&e) 

Pi 

Thermal ion gyroradius(= v t h.i/£li) 

v ) ionization 

velocity averaged ionization cross section 

{<* v ) recombination 

velocity averaged recombination cross section 

MO 

Permeability of free space 
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1 Introduction 


One of the moat interesting and innovative ideas for the production of power in space is to 
generate it by using the geomagnetic field. This can be done with an electrodynamic tether where 
a long cable is swept across the Earth’s magnetic field and current is induced to flow through 
the cable as a result of the motionally generated electric field[l]. In order for the electrodynamic 
tether to be a practical device for generating power in space it is necessary that the load impedance 
and the impedance of the current circuit through the tether ends and the ionosphere be not too 
different. This requires that the tether make efficient electrical contact with the ionosphere. One 
way of doing this is with a metallized balloon which offers a large collection area to the ionosphere. 
This suffers from the disadvantage that the neutral drag on the balloon may be liurge and affect 
the dynamics of the tether. Another means of effecting electrical contact is with the use of plasma 
clouds. Devices which emit plasma clouds for the purpose of facilitating electrical contact in space 
have been named plasma contactors. The expectation is that the plasma cloud will form a large 
interaction region with the ionosphere and allow the motion of charge through the cloud for a small 
potential drop without suffering from the problem of neutral drag associated with a metallized 
balloon. 

The expansion of the plasma cloud and the collection of current has been the subject of both 
experimented and theoretical work. Ground based experiments have proven the basic principle of 
plasma contactors, that they can act as devices to magnify the collection of current from an ambient 
plasma[2]. Theoretical work until recently has focused on the plasma flow near the plasma source 
in the ground based experiments^] and has obtained qualitative agreement. Recently it has been 
realized that the use of plasma contactors in the ionosphere where the cloud can freely expand may 
be different from the ground based experiments and so theoretical analysis has been undertaken for 
space based contactors [4, 5]. These studies show that in order for contactors to work as advertised 
as electron collection devices it is necessary for the plasma cloud to be highly turbulent in order for 
the ionospheric electrons to scatter across the restaining influence of the geomagnetic field. This 
observation is amply confirmed in the ground based experiments. All theoretical studies so far have 
assumed that the plasma flow is dominantly radial away from the plasma source. This assumption 
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simplifies the analysis considerably and enables simple estimates to be made of the current voltage 
characteristic of the plasma cloud[5]. In this paper we relax this assumption and construct a model 
for the anisotropic structure of the plasma cloud across the magnetic field. We show that the 
cloud physics can be divided into three regimes. The first regime is the core region of the plasma 
cloud about the plasma source where the plasma cloud is highly diamagnetic and the electrons are 
effectively unmagnetized either due to exclusion of the geomagnetic field or to turbulent collisions. 
The ions in this region are unmagnetized because the ion gyroradius is expected to be much larger 
than the core region. The core region will be isotropic because there is no preferred direction. The 
second region is where the electrons are magnetized but still experiencing turbulent collisions and 
the ions are unmagnetized. The ion flow here will still be radial. We shall call this the transition 
region. Finally the third region will be where the ions and the electrons are magnetized. This is 
where the plasma cloud will be in contact with the ambient plasma. This region is very important 
to understand because it is here that some of the electrons will be supplied which will eventually 
be collected at the center of the cloud. This region we «hall call the outer shell. This is because 
the major effect of the cloud is to facilitate collection of ionospheric electrons and this outer shell 
provides the collection area akin to the outer surface of a metallized balloon. 

In principle, plasma contactors can be used for both electron and ion collection. However in 
this work we shall concentrate on electron collection. Recent work has dealt with electron emission 
using plasma contactors [6]. The central question we shall address is the use of plasma contactors 
for electron collection in the ionosphere. We are concerned with use of plasma contactors for tethers 
which can draw many amperes of current [1], There are two possible sources for the electrons which 
may be collected by such a device. The first is ionospheric electrons. That is the electrons will be 
pulled in from the far field. The second is by ionization of local neutral gas. In the electric field 
around the collector the newly born electrons will be collected and the ions repelled. We note that 
there are two sources of neutral gets, first the ionospheric neutral oxygen and secondly the neutral 
gas emitted from the contactor itself. It is very important to determine by which mechanism a 
plasma contactor may work since the system implications are quite different for the two possible 
mechanisms. We shall show that if plasma contactors collect most of their electrons from the far 
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field then small plasma contactors are more efficient than large ones which suggests use of small 
multiple tethers for generation of a given current. On the other hand if plasma contactors work 
by ionization of neutral gas then large collectors are more efficient which drives an electrodynamic 
tether system towards long single tethers. We shall address this question of how contactors may 
work in space by constructing a model for the potential in the far field and calculating the electron 
current which can flow in towards the center. We do this in the presence and absence of ionization 
of the neutral gas emitted from the contactor. 

The outline of the paper is as follows: in section 2 we review and further develop the theory 
of the inner core region. In section 3 we consider the outer shell and develop a two dimensional 
model for the potential structure across the magnetic field. In section 4 we solve these models 
numerically and obtain an estimate of the current voltage characteristic of the contactor cloud. In 
section 5 we conclude and suggest the impact of this work on engineering use of plasma contactors. 
We note that there is a large body of work which has been done on electron collection by electron 
beam-emitting vehicles[7]. This work also suggests the importance of ionization to the creation of 
the return current necessary for charge conservation. It has been observed that ampere level return 
currents have been collected for modest vehicle potentials such 100 or 200 Volts. In the case of a 
plasma contactor used for electron collection there are at least two major differences as compared 
with the emission of an electron beam and collection of a return current. The first is that there is 
no perturbing electron beam and hence no possibility of a beam plasma discharge. The second is 
that we are interested in use of a contactor for steady state electron collection and there is evidence 
that it may be possible to collect more current in a transient operation than in steady state[8]. 
These considerations suggest that the physics of electron collection in space by plasma contactor 
clouds may be different from the experience obtained from electron beam experiments. 

2 Review of the core cloud region 

In Figure 1 we show a sketch of what the plasma contactor cloud may look like as used on the 
end of an electrodynamic tether. Near the plasma source the cloud will be very dense and highly 
diamagnetic. In addition the anodic end of the tether will probably be biased to a voltage which is 
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many times the ambient ion energy. These conditions are such that there will be a core region to 
the plasma contactor cloud where the expansion is radial. This can be shown by noting that the 
two directions of anisotropy in the cloud will be the direction of motion of the whole system and the 
direction of the magnetic field. When the plasma is highly diamagnetic or highly collisional then 
the magnetic field is shielded out of the plasma so that this direction of anisotropy is destroyed. 
The other direction of anisotropy is manifested through the motional electric field which is shielded 
whenever the plasma is dense enough so that it can easily polarize. The manifestation, of these two 
directions of anisotropy divides the plasma cloud into three regions. The first region as we have 
discussed is a core region where both the magnetic field and the motional electric field are shielded 
from the plasma. The plasma will expand isotropically in this region and will be determined 
by the boundary conditions imposed by the plasma source. The boundary of this region will be 
where the condition v t /£l e = 1 is satisfied and where the gyrofrequency is defined in terms of the 
actual magnetic field as modified by diamagnetic effects. This condition states that the collisional 
scattering due to either classical or turbulent effects just balances the gyromotion of the electrons. 
This condition should be regared as the most conservative condition for inclusion of magnetic field 
effects. A simple model for this core region was examined in Ref [5]. The second or transition region 
is where the magnetic field is manifested in the plasma but the electric field is still dominated by 
the imposed electric field from the plasma source. In this region the electrons will be magnetized 
while the ions will expand radially under the influence of the imposed potential. The boundary of 
this region will come when the ions have expanded radially a sufficient distance so that they start to 
turn around under the influence of the magnetic field. Hence the boundary will be at approximately 
one ion gyroradius with the gyroradius based on the energy that the ions have gained from falling 
through the self consistent potential. The third region or outer shell will have both the magnetic 
field and the motional electric field manifested. The ions will be magnetized and drift across the 
magnetic field with an & x B/B 2 drift which ultimately is given by the E m x B/B 2 = —V while 
the electrons will flow along the field lines as well as undergo an E x B/B 2 drift. This ensures that 
the ions will ultimately come to rest relative to the ambient plasma. This outer region is where the 
ambient plasma will be contacted. 
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The core region will satisfy the following equations: 
contactor ion density, 


— n n n e ionisation ~ U, M e (<7v) recombination 


neutral density, 


— ~ n n n e (w) ionization + n i n e {& v ) recombination 


potential, 


electron temperature, 


j"; (e(^o - #)) 


m e v e (_J_\ 

en e \4 xr 2 / 


T e dn t ~ dTe 
n e dr dr 


electron density, 



h 

4s - r 2 


E ~ Ei on n n tl e {<Tv) ionigation 


n « — n i H" n amiient • 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where (<rv)i on iMation and (<rw) re combination are given in Ref [5]. In Eq. (4) the electrons are taken to 
be heated ohmically and to lose their energy mainly as a result of ionization. The small loss due 
to electron-ion equilibration is ignored. The plasma sources which are used to produce the plasma 
cloud typically emit cold ions and hot electrons. Hence the ion and neutral temperatures are taken 
to the same with v n given by v n = \/2T,/ m t - and the value of the temperature chosen to be the 
ambient ion temperature. The potential is determined from electron momentum balance( Eq. (3)) 
and quasineutrality is imposed ( Eq. (5)). This means that we are not solving for the potential 
structure in the sheath close to the anode. Rather we obtain the potential drop associated with 
the cloud itself. Hence the potential drops calculated should be regarded as the smallest potential 
drops possible for a given current. 

The ion velocity in the core region will be determined from ion momentum balance. For the 
plasma devices used to produce the plasma clouds the ion are ejected with a velocity which is 
approximately the ion acoustic velocity. In the appendix we show that if the neutral density is 
sufficiently low then the ions will accelerate supersonically away from the contactor source. If the 
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neutral density is high than the ion-neutral friction slows the ions down rapidly and an ion shock 
forms which drops the ion velocity to subsonic values. For parameters of interest the critical neutral 
density is sufficently high(n n » 10 20 meters -3 ) that we shall assume that we are below this neutral 
density everywhere and take the ion to be expanding supersonically. In this case the ion velocity 
will be given by 


( 6 ) 


These equations for the core cloud region must be solved subject to the following 


I = constant 


/ 

while at at r = ro 


n,(r 0 ) = 


4 


n n (r 0 ) = n,(r 0 )( 1 - /,)//,- 


T e = T e (r 0 ) 


and 

dT t T,(r q) 7 e 

dr e 4 jt' 

The electron thermal conductivity is taken to the classical expression[9] 


k = Z.2n t T e /[m t v e ). 


The differential velocity between electron and ions is 

_ I 

VD 4 wr 2 en e ' 

The electron collision frequency is given as 


**• * / «* d" * /,n "I” v acou»iic "I" v Bunt man 


( 7 ) 

( 8 ) 
( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


where for 

v D > c./v^l + v /r e m,/(m e r,)(T e /T ( ) exp(-1.5 - T e /2/T,)) 
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the ion acoustic instability can be triggered and gives[10] 

( 15 ) 

i i Vthe 

and for vd > v t he the Buneman instability can be triggered and gives [11] 

v *acouttic = °-53 (~) 0 61 w P«- ( 16 ) 

The effect of turbulent scattering is included since it has been experimentally observed and theo- 
retically shown to be necessary for efficient operation of plasma contactors[5]. The boundary of this 
core region is taken to be the radius at which the collisionality v t /Cl e = 1 where the magnetic field 
is chosen to be the diamagnetically modified field(that is, B = 0 for /? > 1 and B = B ambient ~ 0 

for 0 < 1). The plasma 0 is defined as 

0 — n e (T« + T,)/(B 2 /2/io). (17) 

The preceeding equations describe the core plasma cloud when the ions are expanding super- 
sonically. The numerical solution of these equations will be discussed in section 4. 

3 Two Dimensional Theory of the Outer Region 

In the outer shell of the plasma contactor cloud where the ambient environment is contacted the 
direction of magnetic field and the direction of motion are the two important directions influencing 
the cloud of plasma which spreads out from a plasma contactor. The magnetic field will influence 
the cloud by allowing it to expand freely along the field but only with difficulty across the field. 
Hence we expect that we can write 


d_ 

dt < dr ± 


(18) 


where l measures along the field § and r±_ measures across B. Eq. (18) is the statement that 
perpendicular variations in all plasma quantities will be much more rapid than parallel variations. 
This statement will only be true where the plasma particles will feel a significant magnetic field as 
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in the outer shell. The plasma cloud projected on the field will be a cigar shape with a possible 
cross-section as in fig. 1. 

The plasma cloud will be described by the state variables n e ,n,-,^,T e and T,- as in the core 
region. Since we expect the plasma elongation we can write for any state variable 

/(fj.,0 - /(rj.) + (/-/> (19) 

where 

Jdl/B I 5 / ( rx,/ ) ( 20 ) 

and 

(/-/)«/ ( 21 ) 

Hence we expect that to lowest order the state variables can be written as independent of the 
position on the field line. The / dtj B average is the correct average that includes the information 
that V • § = 0 everywhere. Physically this is the statement that along the magnetic field, the state 
variables will look the same regardless of location or that variations in physical quantities along 
the magnetic field will be weak compared to variations in perpendicular direction. 

The associated neutral cloud will be expanding spherically and isotropically since the magnetic 
field will not affect it. For temperature of a few eV and for densities of 10 8 cm -s , the mean free 
path is approximately 20 km so that the expansion will be collisionless over the distances of interest 
to us. In that case we can write the neutral density as 

( 22 ) 

where r e is the outer boundary of the transition region of the cloud. 

Since the ion temperature is low on release from most types of contactor so we shall take it as 
constant (and ambient) as in the core region. Since we do not expect extremely large electric fields, 
the plasma will be quasineutral, ft e = ft,- + fi am bitnt • Therefore to complete the state variables we 
need equations for ft<, <j> and T e . These can be obtained by the appropriate averages of the ion 
continuity equation, the equation for current conservation and the equation for conservation of 
electron energy. 
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The / dl/B average can be performed by integrating from the cloud to a distance of the order of 
2(vji/V)r c in each direction. The time 2 r e /V is the transit time of the inner part of the cloud across 
a magnetic field line and v x (the Alfven velocity) is the speed at which information propagates along 
the magnetic field. The distance to integrate to is the maximum distance over which electrons can 
be drawn so that they enter the inner part of the cloud before they sweep past the collector. The 
ion continuity equation is 


dn < ^r> d f n ‘ v ll< 


— + B— (- 
dt dl { 


B 


■) + Vj. • (njVij_) — n g n„ < <TV > ionization fl,' < (TV > recombination (23) 


When the / dl/B average is performed and with the ordering in Eqs. (21) we have 


dfx * 

+ Vx * (ft» Vx t ) = < (TV ionization — Ujft* < W recombination (24) 


In Eq. (24) we have assumed that at the radii of interest to us there is no parallel flow of ions along 
the magnetic field. To complete the ion continuity equation we need the perpendicular ion velocity. 
If the ions are slowly diffusing across the magnetic field then the solution of the steady state ion 
momentum balance[12] gives 

v i± = + k 2 {eS ± + Gi} x | (25) 


where 


and 


_ 1 Ki 

eB 1 + /e* 


k 2 = Kik\ 


(26) 

(27) 


with 


and 



(28) 


Gi = -7\Vln n,- + m,i/, n t; n 


(29) 


The ions are taken to collide mainly with the neutrals (fin)- The first term on the rhs of Eq. (25) 
is the well known Pederson current while the second term is the Hall current. The pickup current 
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has been neglected since in the far field ionization is expected to be small. Since A:,- » 1 (ion mean 
free path long compared to a gyroradius) we can write to lowest order in an expansion in 1/A:,- 

ss + ^V x In n,-) x (30) 

We note that this lowest order expression for the ion perpendicular velocity can only be used in 
the ion continuity equation. For the current balance equation it is well known that the Ex B drift 
for the ions cancels out the electron & x B drift and the higher order terms from Eq. (25) must be 
included to get a nontrivial current density. 

The lowest order perpendicular ion flux arises due to the IxB and diamagnetic drifts. When 
Eq. (30) is substituted in Eq. (24) we obtain 

dft ' T* 

Qt ~ ~ In fii) X B/ B 2 ) = n < av > ionieation — < <TV > recombination (31) 

subject to 

*( r «) = n ie , (32) 

where the density at r = r e can be determined from the solution for the ion density in the two 
inner regions and as |fl| — ♦ oo 

Vjft,- -► 0 (33) 

The potential 4> can be determined from the / dt/B average of the current conservation relation 
(V • J= 0). The average perpendicular electron current is taken to be due to diffusion across the 
magnetic field. Such diffusion can arise due to classical transport (very small) or due to turbulence 
induced transport. We write 

eft 

fi t v tx = — Z? e (V x ft« — - 3 ^-V x ^) + S x B drifts + diamagnetic drifts. (34) 

The electron diffusion coeificent is given by D e = D an + f? c la«*«eoi where the anomalous diffusion 
coefficient is D a n and the classical one is Ddaeiical- The electron flux is composed of an E x B 
drift ( which cancels the ion drift) as well as a diffusion driven flux from either turbulent or classical 
collisions. It is this diffusion driven flux which will give rise to the incoming electron current. In this 
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diffusion driven flux we have neglected a temperature gradient term for the purpose of simplicity. 
The total perpendicular current is 

Jx = «(Mu - n e v e± ) (35) 

where for the ion term we use the full expression given in Eq. (25). 

The field line average of current conservation is 

0 = V -L * L. + f dl/ B^ endpU ^ 

This equation contains a contribution from the parallel electron current flowing through the flux 
tube over which the integration is done. We are looking at the far field and for the magnetic 
field lines which do not intersect the collector. Hence the parallel current from the two endpoints 
cancels out since in steady state there can be no accumlation of current and hence all the electron 
current which is collected from this part of the far field reaches the near field by diffusion across 
the magnetic field. There is a part of the far field which is on magnetic field lines which intersect 
the collector. The contribution from this field aligned part of the far field will be included in the 
next section. 

Since the electron IxB drifts will cancel out with the ion E X B drifts, the Hall terms in the 
ion current are small and from Eqs. (25), (34), and (36) we obtain for <f> 



+ *+* V * — eD e )V ±fl e ) 



+§xv„.Vx («*.!£) = 0 

(37) 

subject to at r = r e 

^ 

(38) 

and in the far field |rj_) — ♦ oo, 

- = 9 x S 

(39) 


In Eq. (37) the first set of terms arises from the diffusion driven mobility terms in the fluxes while 
the second set of terms comes from the density driven diffusion terms in the fluxes. 
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The boundary condition Eq. (39) is the statement that the field far from the cloud is the 
unshielded field seen from a moving frame. This boundary condition on the potential is consistent 
with the condition that in the far field the ions emitted from the contactor will be at rest in a frame 
fixed to the earth. Any other condition on the potential will require some ion motion in the far 
field. In Eq. (37) we have 

A = PW (40) 

The boundary condition 4> c can be obtained as the potential on the outer boundary of the inter- 
mediate zone. 

The electron temperature can be obtained from the averaged energy equation 

3 / q § \ 

2^« ( ~ X -gj ' Vf e ) — — cT ej _ E_i — Ei on h t fi n < (TV >ioni*ation (41) 

where 

r, x = n« v t± (42) 

The electron temperature is determined by a balance between convection, Joule heating and ion- 
ization. This equation is to be solved subject to at r = r e 

T t = T" (43) 

and in the far field |FjJ — » oo, 

V X T« -» 0 (44) 

These equations describe the plasma flow in the far field and are the obvious generalizations of 
the equations for the core cloud. The numerical solution of these equations will be discussed in the 
next section. 

4 Numerical Solutions for the Plasma Cloud 

The equations for the core cloud are a coupled set of first order nonlinear equations in the radial 
variable. These were solved with Gear’s method which is well suited to the solution of such coupled 
equations. The integration was stopped when the collisionality condition was satisfied. In the 
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transition region the ions will still be expanding radially but the electrons will be diffusing across 
the magnetic field. This region is the most difficult to deal with numerically because the two species 
obey different sets of equations. We attempted to model the plasma flow in this region by taking 
the quantities determined by the electrons, namely the potential and the electron temperature to 
be fixed at their values at the edge of the core region and using the radial equations for the ion and 
neutral density throughout the intermediate region. The integration of the ion and neutral density 
equations was terminated when the radial distance was equal to one ion gyroradius as determined 
by the energy of the ions. By this approximation, boundary values for the outer shell of the plasma 
cloud were obtained. This approximation for the intermediate region will give an underestimate 
for the total potential and electron temperature drop and therefore the numerical results reported 
in this paper should be regarded as being optimistic. The important question of the plasma flow 
in this transition regime will be the subject of a future publication. 

The equations for the cross sectional structure of the outer shell are two dimensional partial 
differential equations. Successive Point Over Relaxation(SPOR) was used to solve the elliptic 
current balance equation. The derivatives were approximated by second order finite differences 
and iterations were carried out until the maximum error was less them 10~ 6 between two sucessive 
iterations. Neumann boundary conditions were applied in both the x and y directions. We required 
the potential to match the following boundary conditions in the far field 

V x £ = —V x B 

in the V x B direction and 

V x ^ — 0 

in the direction of motion. The potential at the inner boundary of the region was given by the 
scheme described above for the transition region. 

For the convective equations we used the two-dimensional flux correction method (FCT) of 
Zalesak[13]. The high order scheme was a leapfrog-trapezoidal with fluxes calculated with the 
flux formulae developed by Zalesak[l4] while a donor-cell scheme was used for the lower scheme 
to complete the FCT algorithm. Symmetric boundary conditions were applied on the density 
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distributions and the flux limiter was applied on every iteration. No assumptions were made about 
the symmetry of the solutions so that the computations were carried out through the whole mesh. 

For a plasma contactor the three quantities of interest are the current, potential drop and gain 
as a function of the plasma contactor ion current. The gain is defined as I/I e (ro) and is a measure 
of the current multiplication. The gain can be written in a form which illustrates the two different 
contributions to the total current. From conservation of current we can write 


I = /«(ro) + /«(ro) = Ic(r e ) + U{r e ) 


(45) 


then using Eq. (45) we can write the gain as 

/ 

///c(ro) = 1 + [(ic(r c ) - Ic(r Q ))/I e (ro) + I t (r e )/I e {r 0 )\. (46) 


The first term in brackets on the rhs of Eq. (46) is the gain due to ionization while the second term 
is the gain due to collection of electrons from the far field. 

The electron current flowing inward through the plasma contactor cloud was determined as 
the surface integral of the current density in the outer shell plus the parallel current flowing into 
the core region along the magnetic field. This last term is the contribution from the field aligned 
currents in the far field. Hence we have 


/.= / I.-dS + /|| = /r r + /|| (47) 

Jr=r e 

where iy is the parallel electron current into the core given by Ijj = j\ |. 23 rr| or< and jy, is the electron 
saturation current density associated with the ambient plasma in the far field. In Eq. (47) the 
integration surface for computation of the perpendicular outer current was taken to have a length 
L|| given by the smaller of 

L\\ = 2\j2T t /m t /u t 

and 

L\\ = *r e (v A /V). 


The first length is the twice the mean free path of an electron in the outer region along the magnetic 
field while the second length is twice the length of the Alfven wing generated by the passage of 
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the contactor cloud. This is the rate at which information about the cloud propagates along the 
magnetic field[15]. The perpendicular electron current in the outer shell was obtained from the 
diffusion driven piece of Eq. (34). The electron diffusion coefficent was taken to be the Bohm 
diffusion coeffient(£> a n = ^Bohm = 6.25 x 10 -J T«/B m*/sec where T e is in eV and B is in Telsa). 
The simulations with the Bohm diffusion coefficient model the effect of turbulence in the outer 
shell. Such turbulence will probably arise due to electron runaway leading to Buneman turbulence. 
To model this properly requires a kinetic analysis of the outer shell. We shall defer this for future 
work. We choose Bohm diffusion since this maximises the perpendicular current that can flow 
into the device from the outer shell. A choice of only classical diffusion would give much smaller 
perpendicular currents than we calculate in this section. 

In Tables 1 and 2 we give the total current I against ion current I e for an Argon plasma with 
T e (ro) = 5 eV,/,- = 0.01,0.99 and for a range of I e = 1, 10, 100 Amperes. For / fl <SC 1 no core region 
exists because the central density is too low while for I e > 100 the gas outflow from the contactor 
is high enough to render such devices impractical for long term usage in space. In addition for 
fi = 0.01 and for I e > 100 the central neutral density is so high that ion-neutral friction cannot 
be neglected in the ion momentum equation. For these and all calculations the ambient density 
was n am n ent = 2 x 10 u m~ s and the ion temperature was taken as 0.1 eV. In Table 1 we study an 
Argon plasma contactor cloud which is initially almost totally ionized. The current and potential 
for the outer shell are those that were calculated as the initial conditions for the state equations 
in the last section. The core cloud typically reaches a radius of the order of a few meters and 
therefore allows an parallel electron saturation current of several milliamperes to flow into the 
current collector at the center. We see that the emitted ion current crossing the boundary of the 
outer shell drops slightly with increasing ion current. This is because of electron-ion recombination 
as the electron temperature drops. The electron current from the outer shell drops markedly as the 
ion current increases. This is a reflection of the fact that as the density of the plasma in the outer 
shell increases the density gradient driven piece of the current becomes more effective at balancing 
the inward flowing mobility driven current and therefore the plasma approaches thermodynamic 
equilibrium. The potential drop for the three ion currents is mainly determined by the high density 
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core and is in the range of several tens of volts. We note that for the all the ion currents the electron 
current is being supplied mainly by perpendicular flow across the magnetic field. This was observed 
experimently by Stenzel(8]. From Table 1 we can compute the gain I/I c (ro) and see that small ion 
current contactors have a much larger gain than large ion current contactors. This suggests that 
there is an advantage in terms of efficiency to using several small contactors in series to replace 
one large contactor. In Table 2 we consider an Argon plasma which is initially one percent ionized. 
For small ion currents the system looks like the fully ionized case in Table 1. For large ion currents 
there is significant ionization of the ejected neutral gas that occurs in the body of the plasma cloud. 
This is reflected in the potential drop which increases as the initial ionization decreases. The gain 
for the large ion current contactor case is then much larger than the small ion current case. 

From these results we can conclude that if ionization is ignored then the current that can be 
collected from the far field, even taking it to be driven by turbulence, is small compared to the 
ion contactor current emitted. This is basically because the ambient ionosphere has too low an 
electron density to supply a significant current unaided. By contrast if ionization of contactor gas 
is included then a substantial enhancement of the plasma current is possible. We also note from 
thses results that most of the current gain due to ionization occurs in the high density core. This 
suggests that inclusion of the more complex two and three dimensional effects will give only slight 
improvements on the results from a radial one dimensional model. 

5 Discussion and Conclusions 

We have developed the theory of the physics of plasma contactor clouds used for electron 
collection in space. Such clouds have been shown to consist of three regions, a core region, a 
transition region and an outer shell. Simple equations were developed for the core region and 
outer shell. The plasma structure in the transition regime was modelled approximately. The 
current voltage characteristic of the plasma cloud for a range of ion currents was estimated. 

This work indicates that ionization of neutral gas will be the major source of electrons for the 
current collected to the anode of a plasma contactor. The collection of electrons from the far field 
for ambient conditions will not give substantial increases in the gain of a contactor. 
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This work indicates two major conclusions of significance for engineering use of high current 
contactor clouds. First the biggest current enhancement was obtained with high ion current con- 
tactors emitting a weakly ionized gas cloud. This implies that a single tether with a one high ion 
current contactor will be more efficient than a set of smaller contactors. However most of the gain 
occurs as a result of ionization of the neutral gas emitted from the contactor itself which indicates 
a large mass flow rate for the highest ion current cases. This suggests the second conclusion, that 
if system considerations preclude high mass flow rates and we wish to use a contactor which works 
by collecting electrons from the far field then the smallest ion contactors are the most efficient. This 
suggests that an array of electrodyamic tethers spaced far enough apart so that their contactor 
clouds would not interact would have a much larger gain than one large ion current contactor used 
on a single tether. This can be illustrated by a simple calculation based on Table 1. A single 
electrodyanmic tether using a 10 Ampere contactor would allow an 10.2 Ampere current through 
the tether. An array of ten tethers electrically in series each with a one ampere contactor would 
enable a 12.1 Ampere current to flow though the tethers. An array of multiple tethers is also more 
desirable for use in space for reasons of reliability [1], 

Future work will concentrate on a more careful treatment of the transition zone. This is im- 
portant to understand because if the potential drop in this region is too high then operation of 
electrodynamic tethers will be adversely affected. 
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A Analysis of the Ion Momentum equation 


In a spherical geometry the ion momentum equation with spherical expansion for the neutrals 


IS 


dvi Ti din ni e d<f> , \ { r oY , \ 

where Vin(ro) is the ion-neutral collision frequency evaluated at the center. We shall assume that 
n,- = n e which is a very good assumption since we are concerned with the high density core where 
the ambient ion density can be neglected. If the current driven term in Eq. (3) is small and since 
the numerical calculations show that dT e /dr < T e /n e dn e /dr we write 

din n e 


'fr= T '- 


dr 


(49) 


If we substitute Eqs. (48) into Eq. (49) and take from continuity (ignoring ionization and recombi- 
nation) 

dlnn,- 2 1 dv { . . 

~Tr 7-7,17 (50) 

then we obtain 

i a ., •> 

(51) 


, l.du »75 2 

(U -^^=^ U+ - +4U ' 


i 

2 U " 


l x du _ 

~ v* n r 

where u = v,/c„ r) = Vi n {rd)r/c t> r?o = i/, n (ro)ro/c,, ti n = t > n /c, and we have used the more 
general definition of the ion acoustic velocity as cj = (T e + T,)/m,-. Since from Eq. (50)we have 

— « 7*7, < 52 > 

we require u ~ l/rj* for ij » 1 in order that the density be regular in the far field. It is easy to 
show that the solution to Eq. (51) which satisfies this is for r) » 1, 

“ = jn: 5 ( 53 ) 

*lo + 

where c is an arbitrary constant. Hence we see that in the far field the solution of the full radial 
momentum equation is deretmines by the neutral collisions. An examination of Eq. (51) indicates 
that the equation is singular along the line u = l(i.e., the flow is just sonic) except for the point 
rj = fji, u = 1 where 

jj 1 = fu - “-)• («) 

rjo 1 
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The point 17 = r)i is the point at which the left and right hand sides of Eq. (51) vanish. If Eq. (51) 
is wriiten in the form 


du _ fett+l + feu* _ /(r;,u) 

dr) u-l/u u-l/u 

and if u n < l(which is reasonable since T e /T{ » 1) then f(r), u) = 0 along the line 

»?(«) = y(« ~ «n)- 


(55) 


(56) 


The line r)( u) divides the (u, rj) space into increasing and dcresing trajectories. For rji > r ) 0 the 
u, T) space will look as in Fig 2 where a typical trajectory is sgown. For ions which are emitted 
with velocity just above the ion acoustic velocity the ion-neutral friction slows them down to 
sonic speed. Since the sonic line is singular an ion shock forms which connects the supersonic and 
subsonic trajectories. If the radial distance given by r) 1 is inside the core region then this solution 
is reasonable. For r)\ > r)o the u,r) space will look as in Fig 3. For this case a trajectory which 
starts slightly supersonically can accelerate because the ion neutral friction never gets high enough 
to affect it. The condition on whether the ions can accelerate supersonically is r) 1 < tjq which gives 


Vinjr 0 ) 
2c, 



(57) 


For r 0 = 0.1 meters, an Argon plasma with T e (ro) = 5 eV and T,- = 0.1 eV and an ion-neutral 
cross section of 10 -19 m 2 this gives n„(ro) < 10 2 ° m~ 3 . hence as long as the initial neutral density 
is below this the ions can be taken as accelerating supersonically and there is no possibility of an 
ion shock. 
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Figure Captions 


Figure 

Figure 

Figure 


Possible plasma contactor cloud structure on end of electrodynamic tether 
u versus rj for r)i > t)o 
u versus rj for 171 < rjo 
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/c( r core)(A) 

'll(A) 

r cort (m) 

If*" (A) 

1(A) 

^^core(V) 

A^outer (V) 

0 0 

1 

i 

1.31X10" 3 

0.35 

0.21 

1.21 

16.5 

1.17 


10 

9.99 

8.63 xlO" 3 

0.95 

0.2 

10.2 

27.13 

1.41 

28.54 

100 

99.9 

7 x 10 -2 

4.8 

7.4X10' 2 

100.04 

39.15 

1.16 

40.31 


Table 1: Current and voltage against ion currents, ro=0.1m, T e (ro)— 5 eV, fi— 0.99, 

n ambient = 2 X 10 11 m *, T| - =0.1 eV, Argon 


4(ro)(A) 

/c( r cor«)(A) 

W 

r core( m ) 

I° u<er (A) 

1(A) 

A&or«(V) 

A0 out « r (V) 

0o 

1 

1.008 

1.31 X 10” 3 

0.35 

0.21 

1.22 

16.54 

1.17 

17.71 

10 

10.81 

9.64 X 10 -3 

1.0 

0.16 

10.98 

27.45 

1.26 

28.71 

100 

219 

0.245 

8.98 

4.8 xl 0 -2 

219.3 

48.1 

1.9 

50. 


Table 2: Current and voltage against ion currents, ro=0.1m, T e (ro)= 5 eV, f=0.01, 

nambient = 2 x 10 u m' 3 , T,=0.1 eV, Argon 
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Cloud Shielding Groduolly Breaks 
Down And Cloud Starts Being Swept 
Back Antiparallel To The Direction 
Of Motion 





